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The problem of the dynamo effect for a Kraichnan incompressible helicity-free velocity field is 
considered. Exploiting a quantum formalism first introduced by Kazantsev (A. P. Kazantsev, Sov. 
Phys. JETP 26, 1031-1034 (1968)), we show that a critical magnetic Reynolds number exists for 
the presence of dynamo. The value of the Prandtl number influences the spatial distribution of the 
magnetic field and its growth in time. The magnetic field correlation length is always the largest 
between the diffusive scale and the viscous scale of the flow. In the same way the field growth is 
characterized by a time scale that corresponds to the largest between the diffusive and the viscous 
characteristic time. 
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The study of magnetic fields generated by the turbulent motion of a charged conducting fluid is relevant for several 
astrophysical applications [Q . An inhomogeneous flow of a charged fluid is able to locally produce a magnetic field. 
The advected field in turn generates electric currents in the fluid and these dissipate the magnetic energy because of 
the finite resistivity. Depending on the properties of the flow, magnetic field creation or dissipation can prevail. 
To determine in general terms the presence or not of dynamo is a daunting task. However, as we will see, there are 
particular models which allow for a detailed treatment. 

The evolution of an initially given magnetic field B(r, 0) in an incompressible flow of a conducting fluid is determined 
by the following equations (T| 



where v(r,i) is the velocity field. The magnetic diffusivity k, which is assumed to be uniform and constant, is 
proportional to the inverse of the electric conductivity of the fluid. 

In Eqs. (|l|) the term (v • V)B is a purely advective contribution that preserves the magnetic energy. The stretching 
term (B • V)v acts either as a energy source or as a sink depending on the local properties of the flow. Finally, the 
diffusive term kV 2 B is responsible for the ohmic dissipation at small scales and balances with the inertial terms at 
the diffusive scale r<j. 

The relative importance of the two contributions on the right-hand side of ([l]) is given by the magnetic Reynolds 
number Re m = UL/k, where L denotes the integral scale of the flow and U is the characteristic velocity at such scale. 
Re m can be regarded as a dimensionless measure of the fluid conductivity. For Re m — > the diffusion dominates and 
the magnetic energy density (proportional to B 2 ) always decays to zero in time. In the opposite limit, Re m — > oo, 
the diffusion term is relevant only at very small scales and the magnetic field is almost frozen in the fluid. We can 
expect that at high magnetic Reynolds numbers the flow is able to enhance the magnetic field, producing a consequent 
growth in time of B 2 . The last process is called dynamo effect, referring to the energy transfer from the velocity field 
to the magnetic one. 

The field B acts on the velocity by means of the Lorentz force, which yields a term proportional to (B ■ V)B in 
Navier-Stokes equations. In general, it would be necessary to take into account such feedback action on v. However, 
since we are interested in understanding if the initial generation of the magnetic field is a persistent situation or not, 
we can assume for the initial conditions B 2 <C v 2 and neglect the Lorentz force contribution. Under this hypothesis 
the evolution equations ([!]) are totally uncoupled from Navier-Stokes equations. Following this kinematic approach, 
we thus proceed as if v was an assigned random field: given the initial condition B(r, 0) and appropriate boundary 
conditions, Eqs. (Q) completely determine the magnetic field evolution. 

For the prescribed velocity we refer to the Kraichnan statistical ensemble g], in which v is assumed Gaussian, 
homogeneous, isotropic and 5-correlated in time. The reason for this choice is that analytical results can be obtained 
HQ. A real turbulent flow is characterized by two scales: the integral scale L and the viscous scale 77, at which 
the dissipation term and the transport one balance in Navier-Stokes equations. The velocity structure function is 
supposed to be smooth for r <C rj, to scale as r% (0 < £ < 2) in the inertial range r\ -C r <C L and to approach 
a constant value at scales much larger than L. The parameter £ represents the Holder exponent of the structure 
function and can be thought of as a measure of the field roughness: for £ = 2 the velocity is smooth in space, while 
with £ = we describe a diffusive field. It is well known that the magnetic dynamo can emerge for an helical flow 
due to the a-effect 0JH]. Here we will restrict to a parity invariant statistical ensemble, which does not give rise to 



I. INTRODUCTION 




dtB + (v ■ V)B = (B • V)v + kV 2 B 



(1) 



1 



a-effect. 

Obviously, several generalizations of the model are possible. In Ref. Schekochihin et al. describe the case of a 
d-dimensional smooth velocity field with a generic degree of compressibility. The a-effect is taken into account in 
Refs. |?],|| respectively in the limit of large and small Prandtl numbers. Furthermore, the effect of the Lorenz force is 
considered in Refs. |^,^[ in the limiting case of very large Prandtl numbers. 

The analysis of dynamo effect is made easier by a simple quantum mechanics formulation, first introduced by 
Kazantsev B . Indeed, on account of the (5-correlation in time of the Kraichnan velocity field, the single time correlation 
function for the magnetic field (i?i(x, t)E>i(x + r, t)) can be expressed in terms of a function that satisfies a one- 
dimensional Schrodinger-like equation. The problem of the dynamo effect can thus be mapped into that of studying 
the bound states of a quantum particle in a given potential that only depends on the velocity correlation function. 
In particular, the ground state energy E of such potential will turn out to be the asymptotic magnetic field growth 
rate. The technique used to compute Eq for different quantum potentials is the variation-iteration method described 
in appendix |A|. 

The aim of this paper is to single out the role that the velocity scales L and r\ play in dynamo theory. To this 
purpose it is interesting to study the magnetic field generation as the magnetic Reynolds number Re m and the Prandtl 
number Pr = v j k are varied (y denotes the viscosity of the fluid). Indeed, they are related to the relative importance 
of the characteristic scales in the physical problem by the expressions Re ~ L/rd and Pr ~ (i]/rd)^ ■ 
In Ref. H Kazantsev finally restricts himself to the limiting case of Re m — > oo and Pr — > and proves that dynamo 
can take place only for a velocity scaling exponent in the range 1 < £, < 2. Here we show that in this range of £ the 
characteristic time of the dynamo effect is of order of the diffusive time td — t^/k = 0(|B|/|kV 2 B|) and the magnetic 
field correlation length is of order r^. We also provide a numerical computation of the growth rate vs £. 
Then, we analyze the case of finite Re m and prove that a critical magnetic Reynolds number exists: if Re m is 
sufficiently small, the dynamo does not ever take place, even for 1 < £ < 2. 

Finally, we show that the Prandtl number does not affect the presence of dynamo, but only determines the magnetic 
field correlation length and the characteristic growth time. If Pr < 1, the field B has a correlation length of order 
rd and it grows with a characteristic time scale of order td- On the contrary, if Pr > 1, the correlation length 
is of order 77 and the characteristic time of order t v , where t v represents the viscous time for the velocity field: 
tv = rflv = 0(|v|/|^V 2 v|). 

This paper is organized as follows. After this general introduction to the problem, in section || we define more 
precisely the Kraichnan model and, following Kazantsev ||], we describe the quantum formalism mentioned above. 
In particular we give the Schrodinger equation that is the central point of this quantum approach. In section III we 



revisit the case of infinite magnetic Reynolds number and zero Prandtl number. Then, starting from these results, we 
analyse the effect of finite Re m and we study how a nonzero Pr influences the magnetic field evolution in time. 



II. THE KRAICHNAN-KAZANTSEV MODEL 

In this section we recall in detail the quantum formalism introduced by Kazantsev in Ref. |^] . The random velocity 
field is assumed to be incompressible, Gaussian, homogeneous, isotropic, parity invariant, and (5-correlated in time. 
Under these hypotheses it is completely defined by its correlation matrix 

(v i (x,t)v j (x',t))=5(t-t')V ij (r) 

= 5(t-t')[D ij (0)-S ij (r)} (r = x-x'), U 

where SV, (r) denotes the structure function of the field v. 

The ^-correlation in time of v is an essential property in order to write a closed equation for the magnetic field 
correlation function that, under a suitable transformation, reduces to a Schrodinger-like equation. 

We impose homogeneous and isotropic initial conditions for B. Therefore, on account of the translational and 
rotational invariance of Eqs. (Q), the magnetic field maintains homogeneous and isotropic statistics at every time t. 
Its correlation tensor has thus the form (see, e. g., Ref. p0|) 



(B i (x,t)B J -(x',t)) = Gi(r,t)% + G 2 (r,t)Af. (3) 



Because of the solenoidality condition V • B = 0, the functions G\ and G2 are related by the following differential 
equation 

9Gi 1 d ,„ 2 



dr r 2 dr 



(G 2 r 2 ). (4) 
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The covariance of B is then completely described by a single scalar function, e. g., its trace H(r,t) — 3Gi(r,t) + 
G2{r,t). The dynamo effect will correspond to an unbounded growth in time of H(r,t). 

The correlation function H(r,t) can be transformed into another function &(r, t) that solves the imaginary time 
Schrodinger equation 



~dt 



Id 2 
m(r) or z 



* = (5) 



in which the mass and the potential depend on r only through Su(r). (For the details see the appendix [y and Ref. 
§)• 

To study the dynamo effect it is useful to put in evidence the time dependence of As usual in quantum mechanics, 
we thus expand the 'wave function' '5 in terms of the 'energy' eigenfunctions ^(r, t) — J il>E{r)e~ Et g(E)dE (or 
^K r ?*) — J2E^E(r)e~ Et for discrete energy levels) and get the 'stationary' equation 

1 d 2 ^ 



m(r) dr 2 



+ [E - U{r)]i> E = 0. (6) 



Referring to the meaning of it is clear that an unbounded growth of the magnetic field corresponds to the existence 
of negative energies in Eq. (^). In particular, it is the sign of the ground state energy Eq that determines the presence 
of dynamo and its value eventually represents the asymptotic growth rate of the magnetic field. Indeed, in this case 
it is the ground state ipE e~ Eot that dominates the growth in time. (Recall that the negative energy levels of a 
Schrodinger equation are always discrete). 

By looking at the variational expression for the eigenvalues in (^|) 

E = fmU^ 2 E dr + f(iP' E ) 2 dr 
J rmpgdr 

one can easily conclude that the presence of dynamo effect is equivalent to the existence of bound states for a quantum 
particle of unit (r-independent) mass in the potential V(r) = m(r)U(r) Therefore, in order to state if dynamo 
can take place for a given velocity field, it is sufficient to study the properties of V. 

Having summarized the quantum mechanics formalism for a magnetic field transported by a Kraichnan turbulent 
flow, in the next section we study the dynamo effect for a velocity correlation function that mimics the real physical 
situation. In particular, we numerically compute Eq and describe the properties of the ground state eigenfunction as 
Re m and Pr are varied. From this analysis we are able to obtain information about the critical magnetic Reynolds 
number, the correlation length of the magnetic field, the asymptotic behaviours of its correlation function, and the 
characteristic time-scale of the magnetic field growth. 



III. TURBULENT DYNAMO 

We consider the realistic situation of a structure function Su(r) that scales as r 2 for r -C 77, as expected in the 
viscous range, as (0 < £ < 2) in the inertial range 77 <C r -C L, and tends to a constant value 25^(0) for r ^> L. 
The case £ = corresponds to the diffusive behaviour, while the other limit £ = 2 describes a velocity field that is 
smooth at all scales below the correlation length L. For the other values of £, in the inertial range Su(r) is only an 
Holder continuous function of r with exponent £/2. The parameter £ thus represents a measure of the field roughness. 

An explicit expression for the velocity correlation tensor, which has the desired scaling properties, is, for example, 

2? y (r) = J e ik r Py (k)d 3 k (8) 

with 

^ e -r\h 

- % 2 + £ _ 2)( , +3)/2 ^-(k). (9) 

The solenoidal projector -Py (k) = (<$y — kikj/k 2 ) ensures the incompressibility of the velocity field. 
In what follows we refer to Eq. (S) whenever we show numerical computations that exemplify our conclusions. 
However, it should be noted that our results are general: they depend only on the qualitative properties of Su(r) and 
not on its explicit form. 
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A. Fully developed turbulent dynamo 



We first consider the limiting case of Re m — > oo and Pr — > 0. Under these conditions the diffusive scale T& is in 
the inertial range and the presence of the cutoffs L and i] is neglected: only the scaling behaviour r* (0 < £ < 2) is 
considered for the velocity structure function. 

The general expression of iSy(r) for an homogeneous, isotropic, parity invariant, incompressible field that scales as 

km Sijir) = D ir i \(2 + £)<% - ^] (10) 



n V 
L — >oo 



where the coefficient D\ has the dimensions oilengthP-® Jt ime. 

In this limit the total energy T>a(0) diverges with the infrared cutoff as L^. 

In order to analyze the existence of the dynamo, let us turn to the quantum formulation described above. The 
potential V has the following asymptotic behaviours (see the appendix [§| and Ref. [[| for the complete expression) 



2/r 2 r < r d 

(2-|£-|£ 2 )A 2 r d «r. 



vw~\; n a , 3 ^/_2 _ z:i (n) 



For sufficiently small £ the potential is positive for all r, it does not generate bound states and therefore the dynamo 
cannot take place. For larger £, V is repulsive up to r ~ r d and becomes attractive at infinity (Fig. A quantum 
mechanical analysis based on asymptotic behaviours ( |ll|) allows to establish that £ = 1 is the exact threshold for the 
dynamo effect Q. 

If < £ < 1 the turbulent flow alone is unable to increase the magnetic field and B 2 decays in time. For those values 
of £, the presence of a forcing term in Eq. (|l|) is necessary to obtain a statistically stationary state. See Vergassola 
H for the detailed analysis. 

From now on we restrict to the values 1 < £ < 2, for which the dynamo is present. 
If Eq. (^|) is rewritten in a rescaled form by means of the transformation r — > r/rd, rd — (k/D) 1 ^, it is easy to see 
that the eigenvalues of the energy must take the form 

E = (12) 

where e(£) depends only on the scaling exponent £ and td is the characteristic time of magnetic diffusion. 
We have already noted that the ground state eigenfunction dominates the evolution in time and that Eg is the 
asymptotic magnetic growth rate. We numerically compute eo(£) as a function of £ by the variation- iteration method 
described in the appendix |A[ The quantity eo grows with £ as shown in Fig. [l| When £ tends to one, eo approaches 
zero and the bound states disappear. In the other limit, eo reaches the value 15/2 according to the results of Kazantsev 
An estimation for eo vs £ already appears in Ref. but there the results are limited to the values 1.25 < £ < 2. 
Moreover, the numerical computations in that paper are performed by a variational method based on the particular 
guess r 2 e~@ r for the eigenfunction ipE - This ansatz is correct for r -C r^, but it fails to capture the right behaviour 
for r ^> rd- Indeed, if we insert the asymptotic behaviours ( pd| ) in Eq. (g), we find that, for 1 < £ < 2, ipE(r) shows 
for r ^> r d a stretched exponential decay with characteristic scale r d and stretching exponent (2 — £)/2 (Fig. |3|). The 
variation-iteration method we used (see the appendix |X|) presents the big advantage of not requiring an explicit form 
for ipE - The algorithm provides as results both the eigenvalue and the corresponding eigenfunction. 



From the expressions of ^^(r) we can recover the behaviour of H (r, t) (seethe definition (B2) in appendix |B|) . We 
have that, for r <C rd, H(r,t) is approximately constant, while, if 1 < £ < 2, the magnetic field correlation function 
decays for r > as a stretched exponential with scale rd 

H(r,t)oce-^ r / r ^ 2 ~ i)/2 (rj«r«L). (13) 

The prefactor 



= « (14, 

depends on the growth rate e (£). We can thus conclude that, for 1 < £ < 2, the magnetic field has a spatial 
distribution characterized by structures whose scales are of order r^. 

The cases £ = 2 and £ = 1 have to be treated separately. Indeed, the asymptotic properties cannot be deduced 
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directly from Eq. (^J). 

The smooth case is solved by Chertkov et al. in Ref. jll| by a Lagrangian approach that relates the growth rate to 
the Lyapunov exponents. There is a big difference between the situation of a smooth velocity field and one that is 
just Holder continuous. In the former case the correlation function is found to depend on the spatial coordinate as 
H(r, t) oc r~ 5 / 2 (equivalent to tpE ( r ) ^ rl ^ 2 )> which implies the presence of structures with at least one dimension of 
inertial range size. Actually the magnetic field in the smooth case has been shown to be characterized by strip-like 
objects. 

Finally, the case £ = 1 can be solved exactly. Indeed, the appropriate ground state eigenfunction of Eq. (||) is (recall 
that £ = 1 is the threshold for dynamo and hence Eq = 0) 

Mx) = c VTT- X (-2 X+{ 2 + X )lo g{ l +X )) i {x = r/rdl (15) 
where the constant C is related to the value of H(Q, •) by the relation C = r\ H(0, •). 

If we neglect logarithmic corrections, the asymptotic behaviour of ipo for r 3> r c i is ipo( r ) °< r 1 / 2 , which yields again 
H(r, •) oc r~ 5 / 2 for r » r^. 

The results we have outlined in this section will be useful in the following to describe the general case where the 
velocity energy spectrum has an infrared and an ultraviolet cutoff. Indeed, we will study a structure function that 
for r <C 7] scales as r 2 and so takes the £ = 2 behaviour, while for r ^> L tends to a constant value like in the diffusive 
case £ = 0. 

B. Finite Reynolds effect 

Let us analyze the situation of finite Re m (and zero Prandtl number). The principal fact is that a large scale 
cutoff L appears for velocity field correlations. The diffusive scale r d is again within the inertial range of the velocity 
fluctuations, and the presence of the viscous cutoff can be neglected. The velocity structure function therefore scales 
as for r <C L and tends to T>u(0) for r ^> L. 

Therefore, the potential V behaves as in the previous case for r <C L, while it takes the £, — behaviour for r 3> L 

!2/r 2 r < r d 

(2 - §£ - \e)/r 2 r d «r«L (16) 
2/r 2 L < r. 

The main consequence of a finite Re m is that V is repulsive also at large scales. It is thus clear that, for sufficiently 
high Re m , a potential well is present at scales of order rd- On the contrary, if Re m is too small, the well can be absent 
or anyway not deep enough to generate bound states (see Fig. [J). Therefore, we can conclude that for sufficiently 
small Re m , the dynamo does not take place, even for 1 < £ < 2. 

The effect of a large scale cutoff on the velocity energy spectrum is thus the presence of a critical Reynolds number 
Rem nt \ For Re m smaller than that value the potential V has not bound states or equivalently, on account of our 
quantum mechanic interpretation, the velocity field is unable to favour the magnetic field growth and the ohmic 
dissipation eventually prevails on stretching. 

The dependence of the dimensionless rate-of-growth Eot^ 1 on Re m is shown in Fig. [|for Re m > Rem in the case 

of the scaling £ = 4/3. Notice that, for Re m 3> Re£f lt \ Eq takes the inertial range behaviour Eq ~ eo(£)tg 1 . 

We can again deduce from Eq. (^|) some properties of the function H (r, t). The correlation length of the magnetic 
field is again of order rd and, at r ^> L, H{r, t) shows an exponential decay 

H(r,i)cxe- 7(r/L) (£<r), (17) 

with 7 = E [L 2 /(2k)]- 1 , R = k + V u (0)/6. 

C. Nonzero Prandtl effect 

Finally, we consider the situation of nonzero Prandl number (at infinite Reynolds number). This is equivalent to 
look at the effect of the viscous scale on the dynamo effect. 

If Pr < 1, the diffusive scale r d is in the inertial range, while, if Pr > 1, it lies within the viscous range. The structure 
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function Su(r) scales as r 2 for r <C rj and as for r ^> 77. 

From the previous considerations we can expect for the potential V the same asymptotic behaviours for r — > 00 as in 
the case of Pr = 0. Therefore, the Prandtl number does not influence the presence of dynamo. What is sensitive to 
Pr is the correlation length of the magnetic field, that approximately corresponds to the scale at which the function 
ipE begin its exponential- like decay. When Pr < 1, the potential has nearly the same shape as in the case Pr = 

- \e)/r 2 r d «r 1 *> 

and the correlation length is of order r^. 

On the contrary, when Pr > 1, the potential well is modified by an attractive £ = 2 contribution 

'2/r 2 r<r d 
V{r) ~ { -4/r 2 r d < r < 77 (19) 

,(2-|e-|C 2 )A 2 V«r. 

For these Pr the function ?pE (r) grows as r 2 for r -C rd, as r 1 / 2 in the range rd -C r *C 77 and has a stretched 
exponential decay for 77 <C r. We can thus conclude that, when Pr > 1, the magnetic field correlation length is of 
order 77. 

On account of what we have just seen, we expect that for Pr <C 1 the ground state energy will be proportional 
to the diffusive time: E Q ~ e o • ^ n the other limit, Pr ^> 1, we can predict an approximate expression for 
Eq by a simple scaling argument. Indeed, for large Pr the potential V behaves like in the case £ = 2 and we can 
expect Eq cx D\ (see Ref. [0 for the discussion of the smooth case). Knowing that Sa(r) cx r 2 for r <C 77 and 
Su(r) oc Dir^ in the inertial range, we can match the previous behaviours to obtain D\ oc 77^ 2 . Finally, we recall 
that from dimensional arguments we have 77 ~ (v/Di) 1 ^ . Summarizing the previous considerations, it is easily seen 
that, for Pr ^> 1, the relation Eq cx t^ 1 holds (the time t v is the characteristic one for the velocity diffusion). 
The Prandtl number Pr ~ {t v /td)^^^~ 2 ^ thus influences also the magnetic field rate-of-growth: in presence of dynamo, 
B 2 increases with a characteristic time-scale determined by the largest beetween the viscous and the diffusive time 
(see Fig. |). 

To conclude this section, we discuss a result that emerges from numerical computations: for Pr ~ 1 the magnetic 
growth rate reaches a maximum (Fig. ^). We can easily guess this behaviour, if we refer once more to the Kazantsev 
quantum formalism. For Pr < 1 the £ = 2 behaviour is practically absent in the potential V , while, when Pr 
approaches the value 1, the scale 77 begins to come into play yielding a strongly attractive —4/r 2 contribution at scales 
r d <C r <C 77. The £ = 2 potential is more attractive than that of other £ and the ground state energy increases in 
absolute value. Then, as Pr becomes larger, \Eq\ decreases as explained above. In other words as long as the viscous 
behaviour affects only the potential shape around r^, its only effect is to make the well deeper and so to favour 
the dynamo. When viscousity becomes very large, the level of velocity fluctuations lowers significantly, inducing 
eventually the depletion of the rate-of-growth. 



IV. CONCLUSIONS 



The presence of dynamo effect for a magnetic field advected by a conducting fluid is strongly dependent on the 
properties of the turbulent flow. We have highlighted in the framework of the Kraichnan ensemble the consequences of 
considering a viscous and an integral cutoff for the velocity. Using the quantum formalism introduced by Kazantsev, 
we have found that a critical magnetic Reynolds number exists. By the same method, we have shown that the Prandtl 
number is the parameter that determines the correlation length of the magnetic field and the characteristic time of 
its growth. Finally, we have argued that in the presence of dynamo the magnetic growth rate is maximum fo Prandtl 
number of order unity. As already noted, the previous analysis depends only on the qualitative properties of the 
velocity structure function, so we expect our conclusions to hold for a generic turbulent flow with same statistical 
simmetries and therefore to be relevant for real applications. 
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APPENDIX A: VARIATION-ITERATION METHOD 

For the numerical analysis of Schrodinger equation we make the transformation y = a~ r (a > 1) which maps 
(0,oo) on the finite interval (0,1). (The constant a should be chosen to properly resolve this interval). Eq. (|6j) can 
thus be rewritten in the form 



Cip = XMip (Al) 



where 



C = -(lna) 2 (yjL + ± ) + {U{y) _ Umm) 

\ dy z dyj y 

M = ^ A = E - U mm 

y 



(A2) 



and U m i n denotes the minimum value of U. C and M. are positive-definite self-adjoint operators defining a spectrum 
of eigenvalues A bounded from below and which extends to infinity. Moreover, C is invertible on all functions twice 
differentiable on (0, 1) and vanishing at the boundaries of the interval. Under these hypotheses the varia tion -iteration 
method described in Ref. jl2) provides a valuable tool to compute the lowest eigenvalue Ao of Eq. ( |Al[ ) and the 
corresponding eigenfunction ip . Indeed, let <po be an initial trial function such that i/>o M. <fio dy 7^ and define the 
nth iterate <p n as 

ip„ = C^Mipn-i = (C- l M) n (po. (A3) 

Then, as n is increased, the sequence ip n converges to the eigenfunction ipQ. The nth approximation to Ao is given by 
the following variational expression employing ip n as trial function 

/ tpnCifindy 

A (") = J± . (A4) 

(finMtpndy 



(n) 

The set A form a monotonic sequence of decreasing values, approaching Ao from the above. The advantage of the 
variation-iteration technique is that no expression is required a priori for the function ipQ. We only have to choose 
any guess for initial function (po and then improve the result by iterating the method for sufficiently large n. The 
convergence is more rapid the smaller is the ratio between Ao and the following eigenvalue. 

Finally, for the numerical implementation of the method, we exploited the first order discrete expression of C 
preserving the boundary conditions on tp. If (0, 1) is divided in intervals of length A and yi = iA, we have 



m(y^) fin 1 

c a = -z^( u (w) - c/ ™») + x \ (w-i + 2 Vi + w*+0 if * = 3 ( A5 ) 



{-Vi-i - Hi) if i = 3 + 1 

(Vi-x + 2y< + Vi+x) i£i=j 
{-Vi - Vi+x) if i = j - 1. 



APPENDIX B: THE SCHRODINGER EQUATION IN THE DYNAMO THEORY 

In the present appendix we refer to the notation adopted in the body of the paper. So, the trace of the correlation 
tensor (Sj(x, t)Bj(x + r, t)) will be denoted by H(r, t). 

As a consequence of the velocity <5-correlation in time, H satisfies a closed equation that, under a suitable transfor- 
mation, takes on the form of a one-dimensional Schrodinger-like equation. In order to exploit this fact, let us denote 
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s(r) = Sa(r) and define the following quantities 

s(r) = 
A(r) = «+3(r), 

Then, the function 

V(r,t) = y/n exp (^J 
solves the imaginary time Schrodinger equation 



1 f r s( P ) a . 



Ai(p) 
2M(P) 



Ai(r) = A(r) + 3 K +^. 



^ H{ P ,t)p 2 dp 



(Bl) 



' St" 



i a 2 

^(r) <9r 2 



-ET(r) 



* = 



where 



m 



1 

2A' 
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Ids 1 A 2 
r cir 2r 2 A 



A- 



Ai 

dr V rA 



(B2) 



(B3) 



(B4) 



(See Ref. 0] for the detailed derivation). If we expand \& in terms of the energy eigenfunctions ^(r,t) 
f tpE(r)e~ g(E)dE, we get the stationary equation 



1 d 2 ^ E 
m(r) dr 2 



+ [£? - U(r)]ip E = 0. 



(B5) 



The dynamo effect corresponds to the presence of negative eigenvalues in Eq. ( |B5[) . 

The correl atio n function H(r, •) must tends to a constant value as r — > and decreases to zero as r — ► oo. From the 
definition (jB2| ) we have therefore that Eq. (B5) must be solved with the boundary conditions that ipE(r) vanishes as 
r —* and increases as r — > cxd slowly enough to guarantee that H(r, •) decreases to zero. In particular, if s(r) tends 
to a constant as r — > oo, tPe{t) cannot increase more rapidly than r. 
We consider now the explicit expression 



with 



V ij (k)=D 



-r}k 



(k 2 +i-2)(?+3)/2 y 



r X> tf (k)d 3 k 



P y (k), (a > -1). 



(B6) 



(B7) 



The transverse projector Py (k) = (<5y — kikj/k 2 ) ensures the incompressibility of the velocity field. 
In the limits 77 — 5- and L — > 0, Sij(r) takes the form 

km S i3 -(r) = L>ir« [(2 + - C^? 

L — >oo 



(B8) 



with 



Di 



4™(f )r(-i-o 



C + 3 



-D . 



(B9) 



(The function T is the Euler function) 



If we insert s(r) = 2(£ + 3)Dxr^ in (B4), the transformation (B2) takes on the form 



tf(r,t) = 



H(p,t)p z d P , 



(BIO) 
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while its inverse reads 



H(r,t) 



{2 k- Dx r s (f - 2)) *(r,f)+2r( K + D x r«) *'(r, t) 



2r 2 (K + L»ir«) ! 

For the mass and the potential we obtain the following expressions 

1 



m(r) 



2(K + Dxrty 



(Bll) 



(B12) 



U(r) 



4k 2 + A{t)KD ir S + B{£)D\r 2 Z 
t 2 (k + Dirt) 



(B13) 



with = (8 - 3£ - C 2 ) and 5(0 = (4 - 3£ - §£ 2 ) 
For the sake of completeness we write also the expressions of the trace s(r), which we used to compute Eq respectively 
in the case of finite Re m and in the case of nonzero Pr 



lim s(r) 
?7— >o y 



/ r 2 
VIZ 2 



1 - - 

1 2 

€±i I o 

2 ' 2 ' 



(B14) 



lim s(r) = 8nD ^ ( r(-£) + - ( 1 + ^ ) ' - sin[(l + arctan(-)] I . (B15) 



(The function G denotes the G-Meijer's function of argument r 2 /(4L 2 ). See Ref. |T^] for the exact definition). The 
explicit expressions of the mass and the potential can be derived from (B4). 
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FIGURES CAPTIONS 

Fig. 1: The dependence of the magnetic growth rate eo = ^o^J 1 on the scaling exponent £ in the limit of infinite 
Re m and zero Pr, as computed by the variation-iteration method described in appendix |A|. 

Fig. 2 : The shape of the quantum potential V in the limit Re m — + oo and Pr — > for a value of £ (£ = 2) 
for which dynamo is present and for a value (£ = 0.91) for which there is no dynamo effect. 

Fig. 3: The asymptotic behaviours of the 'stationary wave function' ipE in the limit of infinite Re m and zero 
Pr. The maximum at r ~ r^ determines the magnetic field correlation length. 

Fig. 4: A qualitative picture of the quantum potential shape for Re m respectively above and below the critical 
value Ret U) (K £ < 2). 

Fig. 5: The dependence of the magnetic growth rate on the magnetic Reynolds number for Pr — > and £ = 4/3. The 
numerical computation is performed using expression (|^) for the correlation tensor of the magnetic field. 

Fig. 6: The dependence of the magnetic growth rate on the Prandtl number for £ = 4/3 and in the limit Re m — » oo. 
The numerical computation is performed using expression (P) for the correlation tensor of the magnetic field. 
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FIG. 1. D. Vincenzi 
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FIG. 2. D. Vincenzi 
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FIG. 3. D. Vincenzi 



13 




FIG. 4. D. Vincenzi 
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